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FOREWORD 


This report presents work done by the Defense Logistics Agency Operations 
Research and Economic Analysis Management Support Office (DLA-DORO). This 
study was made possible by the Defense National stockpile Center, Office of 
Planning and Market Research, Messrs. Richard Corder and Franklin Ringquist. 

The resulting product from this study will enable the Defense Logistics Agency 
to more quickly prepare an optimal Annual Materials Plan. Timeliness and 
optimality are essential since the Strategic and Critical Materials 
Stockpiling Act (50 U.S.C. 98 et seq.) requires the President to submit the 
plan to Congress each year. The Executive Branch has delegated preparation of 
the Annual Materials Plan to the Defense Logistics Agency, Defense National 
Stockpile Center. 
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BCBCraiVE SOMffiY 


At the request of the Defense National Stockpile Center (DNSC ), the Defense 
Logistics Agency Operations Research and Econcndc Analysis ttanaganent Seaport 
Office (DLA-DCRO) developed the Annuad Materials Plan Analysis Tool (AMPAT). 
Diis tool quicdcly and acxxnrately builds an Annual Materiads Plan given 
macro-level factors such ais Tninimim and maxiraim acsquisition/di^xssail 
quantities, foreioasted prices, and legislative di^xsead authority. 

Historically, the pr^aration of these plans has beai prinarily manuad, and 
therefore, time oensuming. Ihis becomes irportant vhen one considers the 
.sometimes frequent changes to the plans due to Defaise Department and 
Oongressionad interest and the attendant quick re^xxises required. AMPAT has 
gi>Aen ENSC the c^ability to produce the AMP cadcadations within hours as 
opposed to weeks. 

Initiadly, it was planned to incorporate this tool in a larger oonputer 
program that would aid the ENSC industrial speciadists in developing the 
macro-level factors. AMPAT would be incorporated in the program to lay out 
the optimum levels and timing of acquisitions and disposals given acquisition 
and di^xssad authority and other practical constradnts. i\fter observing an 
actual preparation of an AMP, it became apparent that the DNSC quickly needed 
a tool that solved the problem of oenputing an AMP that satisfied the given 
restrictions. Discussions with the management of ENSC led to an agr eemen t to 
split the development into two parts: first, develop a tool to cxnpute the 
AMP given macro-level factors; and seccmvl, a follcw-on project wcxdd 
enhance AMPAT by vising artificial intelligence concepts to develop the 
macro-level factors. 


Ihis report marks the completion of the first project. Ihe initiad version of 
AMPAT has been developed and delivered to ENSC for immediate use. Ttie report 
recommends that the second phase be started immediately in order to satisfy 
the full requirement. 


lx 


1.1 


BACKQROXniD 


SECTION 1 
INTRODUCTION 


The Annual Materials Plan (AMP) states how much of what materials to add, to 
remove, or upgrade in the Defense National Stockpile (DNS). The AMP also 
states when during the next 3 years such activity should occur. Materials 
found in the DNS are those which would be scarce in the opening years of a 
global conventional war that forces full economic mobilization of the United 
States. Public law creates the AMP which the President submits to Congress. 

The AMP is a relatively short document but is tedious and time-consuming to 
create. It requires coordination with many federal agencies and is affected 
by legislation and politics. It requires understanding the economy and 
political stability of foreign supply sources, international supply and 
demand, and being able to forecast commodity prices. Although these factors 
are volatile and will always be difficult to predict, it is possible to 
mathematically model some important interactions. 

The ability to swiftly build or quickly modify an AMP is important because of 
the scrutiny the AMP receives at the highest levels of the federal government. 
It is a time-sensitive document that often undergoes rapid changes from within 
the Office of the Assistant Secretary of Defense/Production & Logistics 
(OASD/PStL), sometimes with Congressional interest. 

1.2 PURPOSE 

The purpose of this study is to enable the Defense National Stockpile Center 
(DNSC) to swiftly compute or modify an AMP. 

1.3 OBJECTIVES 


The general objectives of this study are to: 


(1) Understand how economics and politics affect an AMP. 

(2) Understand the techniques and philosophies used to build an AMP. 

(3) Build a personal computer based analysis tool that implements these 
techniques and philosophies and which allows DNSC industrial 
specialists to build an AMP. 

(4) Lay the foundation for a later study that uses artificial 
intelligence (AI) concepts to improve the process of developing the 
input parameters to the AMP. 


1.4 


SCOPE 


This study uses gross macro-level factors that affect an AMP and whose 
interactions can be numerically modelled. Also, this study and the resultant 
computer based analysis tool only address the last stage of building an AMP. 
In this stage, the macro-level factors have already been estimated. At the 
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request of the DNSC, the 3 year horizon required for the AMP is extended 
to 10 years for the purpose of the model. 

This study will not address the decision making processes or the derivation of 
the macro-level factors and other inputs to the model. These concerns will be 
addressed in a later project which uses AI and operations research concepts to 
model the decision making process that precedes computing the AMP. 
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SECTION 2 
METHODOLOGY 


2.1 


CONDUCT BACKGROUND INTERVIEWS 


In July 1991, representatives from the ONSC Office of Planning and Market 
Research (DNSC-R), the OLA Operations Research and Economic Analysis Office 
(DLA-LO), and the OLA Operations Research and Economic Analysis Management 
Support Office (DLA-DORO) met to discuss if AI could help ease the task of 
building the AMP. The result of this meeting was a concept paper that said AI 
and operations research concepts could be used. The resulting analysis tool 
would be named the Annual Materials Plan Analysis Tool (AMPAT). 

In December 1991, DNSC industrial specialists attempted a week-long training 
session that would have taught DORO how to build an AMP. Unfortunately, 
urgent taskings from OASO/P&L regarding the AMP completely upset training. 
Although training could not achieve its goal, the rapid-fire taskings from 
OASD/P&L (arising from Congressional actions) dramatically underscored the 
need for the DNSC to rapidly build and modify an AMP. It was decided to split 
the effort into two parts, the first part addressing the most critical need of 
the DNSC: computing the AMP given that certain rates have already been 
decided. The second part would incorporate AI techniques to enhance AMPAT and 
to model the entire decision making process that determines the input factors 
used to compute an AMP. 

2.2 IDENTIFY MACRO-LEVEL FACTORS AND MATHEMATICAL MODEL 

Training did, however, successfully identify the macro-level factors that 
affect how an AMP is built. Furthermore, by viewing how the AMP was computed 
in response to the taskings from OASD/P&L, it became clear that the DNSC was 
unknowingly trying to manually solve a mathematical problem known as a linear 
program (LP). The LP often arises in many practical applications and a 
solution method is well known. 

The data required to solve the LP are global and commodity-specific factors. 
The global factors are the operating and maintenance costs, research and 
development costs, the maximum spending allowed for acquisitions, and the 
upper/lower dollar limits imposed upon disposal sales. Commodity-specific 
factors give the present status of a comnodity. These factors are an item's 
current inventory, how much extra may be currently disposed (i.e., unused 
disposal authority), and aggregate requirements during the planning horizon. 
Commodity-specific factors that vary annually are the upper/lower limits on 
how many units can be bought or sold, the average unit price, and additional 
disposals. 

2.3 IDENTIFY VALUES FOR MACRO-LEVEL FACTORS 

2.3.1 SOURCE DOCUMENTS 

The values that macro-level factors can assume and the manner in which they 
can interact are partly determined by legislative and technical documents. 

The essential character of the AMP is found in the public law that establishes 
the DNS. The Strategic and Critical Materials Stockpiling Act (50 U.S.C 98 et 
seq.) mandates that a stock of strategic and critical materials is to be 
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maintained to decrease and preclude dependence upon foreign sources of supply 
in times of national emergency. The AMP is continually affected by 
legislation and political interests that affect the amount and type of 
materials acquired/disposed. 

Another document that affects the AMP is the Report to the Congress on 
National Defense Stockpile Requirements. This document is submitted by the 
Secretary of Defense to Congress. Recommendations found in this report are 
based on Department of Defense assessments of national needs in conjunction 
with a military conflict scenario developed by the Joint Staff. This scenario 
determines theaters of conflict, military force structure, losses, 
military/industrial and essential civilian requirements, domestic production, 
and supplier country reliability during the first three years of a global 
conventional war of indefinite duration. The requirements used to build the 
AMP arise in the analyses considered in this Report. 

2.3.2 SOURCE MODELS 

There are several sources for determining the values of the macro-level 
factors that ultimately affect the AMP. The Joint Staff baselines military 
requirements by running the Joint Industrial Mobilization Planning Process 
Requirements Module. Requirements are simply the shortfall between supply and 
military demand. 

Civilian and industrial requirements are also baselined to estimate what would 
be available for the production of defense goods. Demands for the civilian 
sector of the economy are computed from data in the National Income and 
Product Account. 

The above three baselined requirements are then input into the Materials 
Defense Economic Impact Modeling System (MDEIMS) to refine requirements. 

MDEIMS consists of several economic models maintained by OASD (Program 
Analysis and Evaluation). Finally, the Stockpile Sizing Module (maintained by 
the Institute for Defense Analyses) is used to report how much of what 
materials should be stockpiled. These results in turn become constraints in 
the linear program that computes the AMP. 

2.4 DEVELOP THE LINEAR PROORAM THAT COMPUTES THE AMP 


Once the input values for the macro-level AMP factors have been established, 
the mathematical model that computes the AMP can b.e solved. Appendix A shows 
the LP that is used in the AMP. Technical readers may obtain a summary of the 
exact LP approach used (Appendix B). The purpose of linear programming is to 
optimize some function (e.g., maximize disposal sales) subject to many 
constraints. The constraints that arise in building an AMP are: 

(1) Each commodity's annual upper/lower limits on what can be bought or 
sold (in units). 

(2) Each commodity's cumulative annual activity (in units) cannot exceed 
the shortfall (for acquisitions) or exceed the surplus (for 
disposals) in requirements. 
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(3) Annually, total dollars spent on acquisitions, operating/ 
maintenance, and R&D cannot exceed the lesser of a spending cap or 
sales from disposals. 

(4) Annual disposals (in units) of a commodity cannot exceed the 
cumulative annual disposal authority for that commodity. 

(5) Total sales from disposals in each year have an upper and lower 
limit. 


2.5 


IDEMTIFY HOW AMPAT WOULD IMPLEMENT LINEAR PROGRAMMING 


Because an LP-solver would be the heart of AMPAT, an April 1992 market survey 
was done to find suitable and affordable LP software. The survey sought: 

(1) Source code that could be tailored for specific roles. 

(2) Object code that could be integrated into other AMPAT software but 
not modified. 

(3) Stand-alone LP solvers that were optimized for speed and efficiency 
but could not be modifed or integrated into other AMPAT software. 

The reason these different areas were considered was that choosing any one of 
these would fix the manner in which AMPAT would solve an LP. The source code 
approach would allow the most flexibility. The stand-alone LP solver would be 
the least flexible but would be faster. Object code stands in the middle; 
although it originates with source code, it has been compiled into machine 
language (and so can't be modified) which can be directly linked with 
AMPAT. 


The survey paper recommended that AMPAT's LP module be created from scratch, 
because commercial source codes could not be found and object code/stand-alone 
solvers were too expensive. Furthermore, DORO had the necessary compilers and 
skill to write their own source code whose object code in turn could be 
integrated into other AMPAT software. 

2.6 PROORAMMINQ DESIGN OF AMPAT 


Although AMPAT would center on an LP-solver, the time consuming phase of 
developing AMPAT was seen to be creating its user.interface. Also, an eye had 
to be kept to the future: it was envisioned that after a later study, AMPAT 
would become a hybrid product that used operations research and AI to address 
the entire process of creating an AMP. It was thought that the future 
development of AMPAT would be largely in the direction of using AI to model 
decision making. 

The three concerns, therefore, that weighed on the initial design of AMPAT 
were its LP-solver, its user-interface, and its future development. To 
address all three required software tools that would collectively allow 
relatively swift progreun development. Choosing such tools, however, can be 
difficult because a different analyst charged with modification and 
maintainence of AMPAT may prefer different tools. 
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In AMPAT's case, the AI language PROLOG and the procedural language C were 
chosen to develop 7VMPAT. The reason for this choice was that the creator of 
AMPAT had several years of experience building user interfaces and 
applications in these languages. Compilers in both languages were also 
available in DORO. The LP-solver would be created in Borland C++ 1.0, the 
user interfaces in PDCPROLOG 3.21, and PROLOG would control the execution of 
AMPAT. Of all the compilers in DORO, PDCPROLOG had the best facilities with 
which to build user interfaces. 

It should be noted that although FORTRAN was an equally attractive language 
within which to develop the LP-solver, no one (to the knowledge of the PROLOG 
Development Center) had ever attempted an integrated PROLOG/FORTRAN product, 
only PROLOG/C. 
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SECTION 3 

THE ANNUAL MATERIALS PLAN ANALYSIS TOOL (AMPAT) 

3.1 GENERAL FUNCTIONALITY 

AMPAT is a personal computer based analysis tool that allows industrial 
specialists of the ONSC to compute different AMPs and to build/maintain the 
associated commodity and general rates datafiles. By creating and modifying 
these databases, the DNSC can swiftly mix and match different rates and 
commodity files to build AMPs for different economic/political scenarios. 

3.2 SPECIFIC FUNCTIONALITY; MAIN MENU 

AMPAT's greetings screen is shown in figure 3-1. From here the user can 
choose to modify or create different types of datafiles, pair a rate and AMP 
file to compute their associated AMP, or obtain an overview of AMPAT. 

Scattered throughout AMPAT are different help texts, thus precluding the need 
for a printed users manual. 

3.3 SPECIFIC FUNCTIONALITY: MASTER COMMODITY LIST 

Since an AMP is built around a set of commodities, AMPAT provides a master 
commodity list from which items are chosen to build an AMP (figure 3-2). Over 
time, the master commodity list will require updating. The user may modify 
the list by changing the name of an item, its unit of measure, current 
inventory, and current disposal authority. Modifications may be done within 
AMPAT. 


Annual Materials Plan Analysis Tool (AMPAT) 


FI 

Help 


F2 

Modify Master 

Commodity List 

F3 

Create/Modify 

an AMP File 

F4 

Global Rates 


F5 

Create Linear 

Program 

F6 

reserved 


F7 

reserved 


F8 

reserved 


F9 

Overview 


esc 

Exit 



Figure 3-1. AMPAT Greeting Screen 
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FORMAT ; Commodity; standard unit; inventory; current disposal authority 

EXAMPLE: bauxite, refractory; short tons; 15232.32; 1900.23 

Aluminum Metal; ST; 57588; 0 

Aluminum Oxide, Abras Grain; ST; 50786; 0 

Aluminum Oxide, Abras Grn. NSG; ST; 118; 0 

Aluminum Oxide, Fused Crude; ST; 249867; 249867 

Antimony; ST; 36004; 36004 

Antimony, NSG; ST; 7; 7 

Asbestos, Amosite; ST; 34005; 34000 

Asbestos, Amosite NSG; ST; 1; 0 

Asbestos, Chrysotile; ST; 9783; 7700 

Asbestos, Chrysotile, NSG; ST; 916; 0 

Asbestos, Crocidolite (DBA); ST; 36; 36 

Baddeleyite; SDT; 15991; 0 

Bauxite, Met. Grade, Jamaica; LDT; 12288824; 12457740 
Bauxite, Met. Grade, Surinam; LDT; 4908512; 5299597 
Bauxite, Refractory; LCT; 276067; 207067 
Beryllium Metal; ST; 353; 0 
Beryl Ore;ST; 16074; 17729 

Beryllium Copper Master Alloy; ST; 7387; 7387 
Bismuth; Lb; 1631406; 300000 
Cadmium; Lb; 6328570; 6328570 
Chromite, Chemical; SDT; 242414; 0 
Chromite, Metallurgical; SDT; 873122; 0 

FIO Save and Return to MAIN ESC Abort 


Figure 3-2. Coosodity Haster List 

3.4 SPECIFIC FUNCTIONALITY; CREATE/MODIFY AN AMP FILE 

An AMP file is simply a collection of commodities (chosen from the master 
commodity list) whose individual data values have been established in 
sufficient detail to build an AMP. In this AHPAT module, users may retrieve 
an existing AMP file to add/delete commodities, modify their values, or build 
a new AMP file from scratch. Figure 3-3 shows both the menu from which an 
existing AMP file can be chosen and the prompt to create a new file. 

- Hit escape to abort - 

93B.AHP 93C.AMP 

93E.AMP NU93.AMP 

- hit escape to abort - 

Name of new AMP File » 


93A.AMP 

93D.AMP 


Figure 3-3. JUiP File Selection Menu 

Shown is rejection of current choices in favor of building a new AMP file 
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After a file has been specified, another menu appears from which the user 
chooses the commodity to be modified or added (figure 3-4). Commodities that 
already appear in the AMP file are marked with an asterisk. After the 
commodity is chosen, its data values screen appears (figure 3-5). The user 
simply enters or modifies these values which appear: inventory, requirement, 
annual activity levels/prices, existing disposal authority, current disposals, 
and the annual requested disposal authorities. The user may also attach 
narrative remarks to the commodity which explain the established values or 
which may prove helpful to a decision maker. The user may specify up to 10 
years worth of data. 

- ESC saves/goto MAIN; D:\93A.AMP - 

♦Aluminum Metal 
♦Aluminum Oxide, Abras Grain 
♦Aluminum Oxide, Abras Grn. NSG 
Aluminum Oxide, Fused Crude 
Antimony 
Antimony, NSG 
♦Asbestos, Amosite 
♦Asbestos, Amosite NSG 
♦Asbestos, Chrysotile 
♦Asbestos, Chrysotile, NSG 
♦Asbestos, Crocidolite (DBA) 

Baddeleyite 

♦Bauxite, Met. Grade, Jamaica 
♦Bauxite, Met. Grade, Surinam 
Bauxite, Refractory 
Beryllium Metal 
Beryl Ore 

Beryllium Copper Master Alloy 
♦Bismuth 
♦Cadmium 

♦Chromite, Chemical 
♦Chromite, Metallurgical 
♦Chromite, Metallurgical NSG 


Figure 3—4. Listing of Commodities in Chosen AMP Pile 
(♦ indicates items currently in file) 
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Commodity »Bauxite, Met. Grade, Jamaica 
Inventory >>12288824 
Requirement »0 

(A)cquisition or (0)ispo8al »D 


DNSC Requested 


— ACTIVITY 

Min 

LEVEL — 

MAX 

Unit 

Price 

Disposal 

Authority 

500000.0 

600000.0 

20.00 

0.0 I 

500000.0 

600000.0 

20.00 

0.0 1 

500000.0 

600000.0 

20.00 

o 

• 

o 

500000.0 

600000.0 

20.00 

0.0 

500000.0 

600000.0 

20.00 

0.0 

500000.0 

600000.0 

20.00 

0.0 

500000.0 

600000.0 

20.00 

o 

• 

o 


Existing Current 

Disposal Year 

Authority Disposals 

12457740.0 0.0 


FI Help F2 Store Item F3 Remarks F4 Deleteltem ESC Quit 


Figaro 3-5. Commodity Entry Screen 


ECIFIC FONCTIOKALI 


I I I M H 


RATES SC REEN 


An AMP file containing commodity information is only the first of two 
components needed to compute an AMP. The second component is a datafile 
containing annual fiscal rates which affect the overall activity of all 
commodities. Figure 3-6 shows both the menu from which an existing global 
rates file can be chosen and the prompt to create a new file. 

Figure 3-7 shows the screen which displays the rates stored in a retrieved 
rates file and from which the user may modify values. The annual values are 
the spending cap, operating and maintainence budget, research and development 
budget, and the minimum/maximum on disposal sales. As before, the user may 
attach narrative remarks to a global rates file to guide future decision 
making. 
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93A.RAT 

93D.RAT 


- hit escape to abort 

93B.BAT 93C.RAT 

93E.RAT MAC.RAT 


- hit escape to abort 

Name of new Global Rates File » 


Figure 3-6. Globel Rates Selection Menu 
Shown is rejection of current choices in favor of building a new rates file. 



Global Rates 

File: 

D:\MAC.RAT 


(in millions) 





$Spend 

$ Spend 

$Spend 

$Sales 

$Sales 



Cap 

O/M 

RfiD 

Floor 

Cap 

1993 


150.00 

32.00 

25.00 

150.00 

300.00 

1994 

1995 

1996 

1997 

1998 

1999 

2000 

2001 

2002 


600.00 

32.00 

25.00 

80.00 

600.00 

1 

FI Help F2 

Store 

F3 Remarks F4 

DeleteFile ESC Quit 



Figure 3-7. Global Rates Entry Screen 
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3.6 


SPECIFIC FUNCTIONALITY: COMPUTING AN AMP 


After an AMP and global rates files have been chosen, the stage is set for 
computing an AMP. This module simply combines commodity-specific information 
and global rates to create a linear program which is then solved. The solution 
of the LP in turn becomes the levels of acquisition and disposal reported in 
the AMP. These levels are computed for each commodity and for each year. 

Figure 3-8 shows the screen from which global rates and AMP files are chosen 
from. If both files begin with the same year, the user is then asked to 
specify the year which closes the planning horizon. If there are global rates 
for each year in the planning horizon and if each commodity has data for each 
of these years, then an AMP can be built. Otherwise, the user is told either 
that some years lack global rates or that some commodities lack data. The 
user will then return to the module in which the affected global rates or AMP 
file can be modified and made whole. 

Otherwise, a summary is displayed for the user to review (figure 3-9). If the 
user accepts this data, then the underlying linear program will be built and 
solved. If the data is rejected, then the user may modify the appropriate 
global rates and/or AMP files and return again to compute an AMP. The 
solution to the linear program is displayed (figure 3-10) in a table that 
gives the annual levels of activity and associated dollar value for each 
commodity. AMPAT will also identify if other optimal solutions exist but will 
not identify them. 

There is one last thing to note. Even though the user may validate the 
commodity and global rates files, it still may not be possible to build an AMP 
due to some contradiction in this input data. For example, if a disposal 
item's cumulative minimum activity exceeds its surplus, then the user must 
modify the minimum activity of the commodity. AMPAT will automatically check 
the accepted data and alert the user to any one of several contradictions 
found (Appendix A). 
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93A.RAT 

93D.RAT 


choose global rates file - 
93B.RAT 93C.RAT 

93E.RAT MAC.RAT* 

- choose AMP file 


93A.AMP 

93D.AMP 


93B.AMP 

93E.AMP 


93C.AMP 

NU93.AMP* 


Start Year is 1993 


Stop year »1994 


Figure 3-8. Building a Linear Program: Pairing Global Kates/AMP Files 

and Specifying Time Period 
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Figure 3-10. Results of Solving Linear Program: Levels of Commodity Activity and Associated Costs to Report in AMP (partial listing) 













3.7 BENEFIT OF USING AMPAT 

AMPAT empowers the DNSC to swiftly answer the frequent and urgent taskings 
from OASD/P&L for a new AMP. By retrieving data from a user-created library, 
the user of AMPAT may quickly establish the characteristics of an 
economic/political scenario for which an AMP will be computed. Once the 
scenario has been set, its underlying mathematical model is quickly built and 
solved. AMPAT takes the solution and builds the AMP which (after review) may 
lead to other scenarios and AMPs. The feedback loop between OASD/P&L and the 
DNSC is thus greatly accelerated and may lead to a better "final" AMP being 
submitted to the Executive Branch. 
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SECTION 4 
CONCLUSION 


This study reached the following conclusions: 

(1) A mathematically based approach to building an AMP does exist. 

(2) Such an approach would use well-known solution methods that would 
quickly determine the optimal levels to report in an AMP. 

(3) Errors that may otherwise arise from using multiple methods to 
compute an AMP are reduced by using a single analysis tool such as 
AMPAT. 
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SECTION 5 
RECOMMENDATION 

This study reconunends that AMPAT be: 

(1) Used to compute the Annual Materials Plan. This will allow the 
Defense National Stockpile Center to submit more timely AMPs as 
mandated in public law. 

(2) Incorporated in a follow-on system that uses AI and operations 
research concepts to model the decision making process that leads to 
computing an AMP. 
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APPENDIX A 

LINEAR PROGRAMMING FORMULATION FOR THE ANNUAL MATERIALS PLAN 


A-1.1 DEFINITIONS 


This appendix establishes the mathematical model—a linear program—whose 
solution is equivalent to computing an Annual Materials Plan. The model's 
commodity-specific factors for commodity i in Year j are: 

- a^^j requested disposal authority 

aj^O current remaining disposal authority 

equals current disposal authority less disposals to-date 

- L^j lower activity bound 

Pij unit price 

upper activity bound 
x^j activity level 

shortfall or surplus, i.e., requirement - inventory 
The model's global rates for Year j are; 

Cj sales cap 

- Fj sales floor 

Sj spending cap 

Oj combined O&M/R&D budget 

Lastly, define these index sets: 

- A acquisition commodities 

- D disposal commodities 

- Y planning horizon [Yearl, YearN] 

A-1.2 STATEMENT OF THE ORIGINAL LINEAR PROGRAM 


The initial objective function of this linear program is: 

maximize £ ^ij^ij 
i e A U D 
jeY 

where the cost coefficients Cj^j are specific to the type of LP solved. The 
initial constraints of the LP are: 

a. Xj^j e for each i e A U D and jeY (activity bounds) 

b. £ Xj^j s |^il> fot* each i e A U D (shortfall/surplus bounds) 

jeY 
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Oj + E Pij^ij S E Pij^ij' each j c Y 

ieA ieD 


(spending < sales] 


d. E Xj^j < E a^j, for each i c D and k e Y (disposal authority limits) 
j<k j<k 


e. Oj + E PijX^j < Sj, for each j e Y 
ieA 


(acquisition spending cap) 


f. E e [FjrCjl, for each j c Y 

ieD 


(sales floor and cap) 


STATEMENT OF THE REVISED LINEAR PROGRAM 


To obtain the LP used in AMPAT's LP-solver, transform the above LP so that 
all variables have a lower bound of zero. This transformation will yield the 
initial basic feasible solution needed to begin a Phase I optimization 
(Appendix B, Section 1.5) and is accomplished by the substitution: 

*ij ~ ^ij Yij 

where y^j e [0,Uj^j-L£j ]. The resulting LP is thus; 

maximize E c^jy^j 
i € A U D 
j€Y 


subject to: 


yj^j e (0,U£j-Lij], for each i « A U D and j e Y 


b. E y^j < 15^1 - E , for each i e A U D 
jeY jeY 


^ PijYij - ^ PijYij ^ -Oj + E PijLij - E PijLij, for each j e Y 
icA icD ieD icA 


d. E y^j < E a^j - E , for each ieD and k c Y 
j<k jsk j<k 

e. E PijYij s Sj - Oj - E PijLij, for each j c Y 

icA ieA 


f. E PijYij s Cj - E PijLij, for each jeY 
ieD ieD 


g. -E PijYij s -Fj + E PijLij, for each j c Y 
ieD ieD 
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A-1.4 


SIMPLE INSPECTION FOR INFEASIBILITY 


It is possible that a linear program built to compute the AMP may not have a 
solution, i.e., it is infeasible. Infeasibility may arise for many different 
reasons, e.g., specifying a minimum level of acquisition that exceeds the 
spending cap. Although it is generally difficult to identify an infeasible 
linear program a priori, AMPAT does perform several rudimentary checks for 
infeasibility and alerts the user when any of these are encountered. 

To begin with, AMPAT only accepts positive real inputs. Thus all the 
variables in A-1.1 (except for 6j^) will be positive. Further inspection of 
item data and global rates are made to ensure that: 

a^o - 0/ i e D (current disposals valid) 

- Lj^j < U^j, ieAUD, j«Y (activity levels are valid) 

Sj, Fj < Cj, j e Y (sales cap exceeds spending cap and sales 

floor) 

- Oj < min(Sj,Fj), j e Y (O&M and R&D costs exceed neither the 

spending cap nor sales floor) 

- ^ ®il ®i0' i e D (first year minimum disposal level doesn’t 

exceed first year cumulative disposal 
authority) 

After these item-level inspections, AMPAT next combines factors across the 
planning horizon to search for infeasibility. These checks ensure that: 

- £ U^jPij s Fj/ j € Y (maximum disposal sales exceed sales floor) 
icD 

- Ej Lij s |5i|, i e A U D (an item's minimum acquisition/disposal 

across the planning horizon doesn't 
exceed its shortfall/surplus) 

- £ L^jPij £ Cj, j e Y (minimum disposal sales in a year does 

icD not exceed the sales cap) 

- £ LijPij — min {Sj , £ Uj^jp^j} — Oj/ j ^ ^ 

isA icD 

(total expense for minimum acquisitions in a given year does not 
exceed what remains after subtracting o&M/R&D costs from the smaller 
of the annual spending cap and maximum possible disposal sales) 
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Finally, each disposal item is checked to ensure that its requested disposal 
authorities will support its minimum activity levels: 


let RemOAj^j be the remaining maximum disposal authority for each 
disposal item i £ D at the end of Year j £ Y 

— RemDA^j = RemDAi^j_l + a^j - L^j 
RemOAj^Q = aj^o 

- if any RemDAij < 0, then the LP is infeasible 
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APPENDIX B 

SIMPLEX ALGORITHM USED IN AMPAT 

B-1.1 INTRODUCTION 


This appendix gives the mathematical background of the simplex algorithm used 
in the AMPAT LP-solver. In a nutshell, this version of the simplex uses 
known variable upper bounds to reduce the number of constraints considered, 
and uses a single artificial variable in a Phase I optimization to find an 
initial basic feasible solution to the given LP. It is assumed that the 
reader is familiar with elementary linear algebra/matrix theory and has been 
exposed to the vocabulary surrounding the simplex algorithm. 


B-1.2 CANONICAL FORM 


Before an LP can be solved by AMPAT, it must first be pre-processed into a 
required form before solution is attempted. The LP first seen by AMPAT must 
have the following canonical form: 

maximize c'X [1] 

subject to 
AX S b 
0 S X S U 

where c, b, X, and U are columnar vectors, ’ is the transpose operator, and A 
is an mXn matrix of real numbers. 


B-1.3 STANDARD FORM 


Next, AMPAT converts the above canonical-LP into the standard form: 

maximize c'X [2] 

subject to 
AX = b 
0 s X S U 

through the introduction of positive slack variables into the previous 
inequalities. Note that although the same variable names are used, the 
matrices in the standard and canonical LPs are not the same. At this point, 
only the form of the LPs are being established. Also note that the slack 
variables have no upper bound, unlike the variables in [1]. In this case, 
the appropriate entry in U for a slack variable is m. 


B-1.4 USE ONE ARTinCIAL VARUBLE 


The next step is to transform the above equalities so that all of the right 
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hand sides are positive, the coefficients of the slack variables remain equal 
to 1, and that only one artificial variable has been introduced. This is 
accomplished by the following: 

let AX = (V I]*X = b, where I is an identity matrix 

- let (V" l”]*x = b" be the cross-section of AX = b all of whose right 
hand side values are negative 

let S be the index set of rows in this cross-section 



Ej + Sj = bj is replaced by Ej - Ej^ + Sj - sj^ = bj - bj^ 

Notice that bj - b]^ is positive because b]^ < bj by design. Also note 
that the coefficient of the slack variable Sj remains equal to 1. 

finally, replace Ej^ + sj^ = bj^ by -Ej^ - sj^ + a = -bj^, where a is the 
sole artificial variable. 


B-1.5 FINAL FORM OF LINEAR PROGRAM 


After the above steps have been taken, the standard-LP of equation [2] has 
been replaced by an LP of the form: 

maximize c'X [3] 

subject to 
AX = b a 0 
0 S X s U 

where A is an mXn matrix and by construction, has m columns that comprise a 
standard bases for Euclidean m-space, and where a member of this bases 
corresponds to the artificial variable. Note that an initial basic feasible 
solution exists to the above LP: values of basic variables correspond to 
entries in b and non-basic variables equal zero. It is to obtain this 
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initial basic feasible solution that the transformation found in Appendix A 
Section 1.3, was done. The pre-processing phase is now complete and the 
above LP is ready to be solved by AMPAT. 


B-1.6 PHASE I OPTIMIZATION 


If indeed an artificial variable exists in the above LP (as will usually be 
the case), then it is necessary to drive its value to zero to obtain an 
initial basic feasible solution to the original problem given to AMPAT. To 
do this using Phase I optimization, simply change the above objective 
function to: 


maximize -a 

and solve the new LP using the procedure outlined below. If the optimal 
objective value is nonzero, then a cannot be driven to zero and the LP given 
to AMPAT doesn't have a solution. Otherwise, if the artificial variable is 
non-basic, then proceed with Phase II optimization. However, if the 
artificial variable is basic, try to swap it with a non-basic variable and 
then proceed with Phase II. If this isn't possible, then simply recognize 
that redundancy has been discovered and proceed with Phase II. In all cases, 
before going to Phase II optimization, drop the artificial variable from the 
LP and do not ever consider it again. 


B-1.7 THE SIMPLEX ALGORITHM 

B-1.7.1 INITIAL LINEAR PROGRAM 


Because both Phase I and Phase II optimizations follow the same procedure, we 
drop the distinction between the two and discuss how to find the optimal 
solution of a general LP using the simplex algorithm. Let the general LP be: 

maximize c'gXg c'gXg 

subject to 


AgXg + AgXg = b [4] 

0 ^ *B,R ^ ^ 

where B, R are index sets to the basic and nonbasic variables respectively. 
Since Ag is invertible by definition of B, multiply equation [4] by the 
inverse of Ag to obtain: 


Xg + Ag-lAgXg = Ag'lb 


[5] 


B-7 




B-1.7.2 


TRANSFORM LINEAR PROGRAM/CHOOSE ENTERING NONBASIC 


Solve for X 3 in the above equation and substitute into the above LP to obtain 
the equivalent LP: 


maximize [ 0*3 - c' 5 Ab"^A 3 ] * Xr [5.1] 

subject to 

Xb + Ab ^AjjXb = AB~^b 
0 - ^B,R - U 

Now choose the nonbasic whose per unit change increases the objective most 
(cost is synonymous with objective coefficient): 

of all nonbasics at lower bound, let t have the largest positive cost 

of all nonbasics at upper bound, let 12 have the least negative cost 

let ENTER be the nonbasic from {^,(2} whose cost has the largest 
absolute value 


if ENTER is Tull, then optimality has been reached 


Now ask the question "by how much can ENTER change such that Xb,r remains 
feasible?” Lot this max allowable change be 6 . From [5]: 

* AB"^b - Ab'^ArXr 

= AB"^b - y*XR [6) 


Let 


XlR = Xr 1- D 


[ 6 . 1 ] 


where 5 is the only nonzero in D and appears in the position that corresponds 
to ENTER. From [ 6 ], the new Xr is: 

XIr = AB“^b - y*xiR 

= Ao'^b - Y*(Xr + D) 

= A^-lb - X-Xp - YD 

= ' 2 ~ ®*XeNTER [ 1 


where Yr^j^Br is the column in the matrix Y that corresponds to the variable 
ENTER. Note that YD = 5 *Yehter only because 6 is the sole nonzero in D and 
appears in the row that corresponds to the variable ENTER. 


B-1.7.3 CHOOSE EXITING BASIC 

It remains to compute 5 and to determine which basic variable gets replaced 
by ENTER. Per the original LP, 0 < Xr^r s U and so: 

0 < XIr S Ur 

0 S XIr < Ur [7.1] 
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Thus from [7); 


0 S Xb - 5*yENTER - 


which implies: 


~ - **^ENTER - *B 


Let Venter 


[Vi], Xb = 
6 < min 


(Xj^J, and Ub = From [8]: 

{ min { Xi / Vi s Vi > 0}, 

min { [x^-ui] / Vi : Vi < 0} 

} = al (pushes a basic to its extreme) 


S > max { max { 
max { 
} = o2 


Xi / Vi ' Vi < 0>' 

[Xi-Ui) / Vi : Vi > 0} 

(pushes a basic to its extreme) 


[8] 


[9] 


[ 10 ] 


But since only one nonbasic variable changes (namely, xenter) that this 
variable has an upper bound: 

0 < Xenter + 5 ^ 'ienter 


from which it follows: 


-Xenter ^ ^ ^ Center ~ Center 
But since nonbasics are always at their extreme values: 

0 £ 6 £ Center (enter at lower bound) [11] 

““ENTER 5 5 £ 0 (ENTER at upper bound) [12] 

Thus, C9-C12 imply: 

6 = min(ol,UENTER ) ENTER at lower bound 
5 = max{a2,-UENTER)^ ENTER at upper bound 

Although we have an entering nonbasic ENTER, it remains to determine if it 
will indeed replace an existing basic. If 6 arises from [9]-[10], then the 
leaving basic variable is the x^ whose value determines 5. The optimization 
process resumes anew starting at equation [4] (revised in [5.1]) with new 
sets B and R, or if B and R remain unchanged, then at equation [6] with the 
new values for the nonbasic variables (given in [6.1]). 
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